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Spin-transfer ferromagnetic resonance (ST-FMR) in symmetric magnetic tunnel junctions (MTJs) 
with a varied thickness of the MgO tunnel barrier (0.75 nm < tM S o < 1-05 nm) is studied using the 
spin-torque diode effect. The application of an RF current into nanosized MTJs generates a DC 
mixing voltage across the device when the frequency is in resonance with the resistance oscillations 
arising from the spin transfer torque. Magnetization precession in the free and reference layers of the 
MTJs is analyzed by comparing ST-FMR signals with macrospin and micromagnetic simulations. 
From ST-FMR spectra at different DC bias voltage, the in-plane and perpendicular torkances are 
derived. The experiments and free-electron model calculations show that the absolute torque values 
are independent of tunnel barrier thickness. The influence of coupling between the free and reference 
layer of the MTJs on the ST-FMR signals and the derived torkances are discussed. 



I. INTRODUCTION 



High density magnetic random access memories can 
be implemented using current-induced magnetization 
switching (CIMS) Q which is caused by interactions be- 
tween spin-polarized current and the magnetization of 
the free layer (FL) in magnetic tunnel junction (MTJ) 
cells. This phenomenon is called the spin-transfer-torque 
(STT) effect @, ij. Moreover, STT is utilized in MTJ 
nano-oscillators that generate signals in the GHz fre- 
quency range [4-6]. In order to optimize MTJ param- 
eters, so that they can compete with existing memory 
and microwave technologies, it is necessary to fully un- 
derstand STT. The spin-torque diode effect enables quan- 
titative measurements of STT parameters @~H1- I n this 
work, we use the spin-torque diode effect to investi- 
gate the dependence of in-plane and perpendicular spin 
torques on MgO tunnel barrier thickness. The tunnel 
barrier determines the transport properties of the device, 
as it affects the tunneling magnetoresistance (TMR) ra- 
tio, the resistance area (RA) product and the coupling 
between the FL and the reference layer (RL). We show 
that the spin-torque ferromagnetic resonance (ST-FMR) 
spectra contain a double resonance mode for very thin 
MgO barriers due to strong ferromagnetic interlayer cou- 
pling. Moreover, the in-plane and perpendicular spin- 
torques do not depend on MgO barrier thickness, in 
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agreement with free electron models [10j . 



II. EXPERIMENTAL 

The MTJ stack with a MgO wedge tunnel barrier 
was deposited in a Singulus Timaris cluster tool sys- 
tem. The multilayer structure consisted of the follow- 
ing materials (thickness in nm): Ta(5) / CuN(50) / 
Ta(3) / CuN(50) / Ta(3) / PtMn (16) / Co 70 Fe 30 (2) 
/ Ru(0.9) / Co4oFe 4 oB 2 o(2.3) / wedge MgO(0.7 - 1.1) 
/ Co 4 oFe4oB 2 o(2.3) / Ta(10) / CuN(30) / Ru(7). The 
slope of the MgO wedge barrier was approximately 0.017 
nm/cm. The deposition process was similar to the one 
used in our previous studies [TJ Ell ■ After thin- film de- 
position, three different parts of the sample were selected 
for patterning into nanometer size pillars (later in the pa- 
per referred to as SI, S2 and S3, see Table Q]for details). 
Using a three-steps electron beam lithography process, 
which included ion beam milling, lift-off and oxide and 
conducting layers deposition steps, nanopillars with an 
elliptical cross-section of 250 x 150 nm were fabricated. 
The pillars were etched to the PtMn layer. The elec- 
tric leads to each MTJ nanopillar consisted of coplanar 
waveguides which were designed to match an impedance 
of 50 Ohms. To ensure good RF performance, the over- 
lap between the top and bottom leads was about 4 /im 2 , 
which resulted in a capacitance of less than 1 x 10~ 14 
F. Each set of MTJs with a constant MgO tunnel barrier 
consisted of 10 - 15 nanopillars. 

ST-FMR measurement were conducted in a frequency 
range from 2 to 12 GHz. In these experiments, the appli- 
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TABLE I: Summary of static parameters of the prepared MTJ 
nanopillars. 

Sample No. MgO thickness TMR RA product Hs 





(nm) 


(%) 


{Qfim 2 ) 


(Oe) 


SI 


1.01 


170 


9.6 


-21.7 


S2 


0.95 


165 


6.24 


-3.7 


S3 


0.76 


110 


2.86 


47 



cation of an RF current to an MTJ generated a DC volt- 
age (also called mixing voltage V m i x ) across the device, 
when the current frequency was brought into resonance 
with the resistance oscillations arising from the STT. The 
MTJs were placed in an in-plane magnetic field at an an- 
gle of f3 = 70° with respect to the easy magnetization 
axis (except for the case presented in Fig. IHfb)), so 
that a large variety of angles 9 between the junction's 
FL and RL could be obtained. We estimated 8 from the 
assumption, that the resistance R of the MTJ changes as 
follows: 



cos(0) 



Rap + Rp 



R 



Rap - Rp 



(1) 



where Rap and Rp are the resistance of the MTJ for 
an antiparallel and parallel alignment of the FL and RL 
magnetization, respectively. In order to obtain the clear- 
est STT results [l3[ , the strength and angle of the exter- 
nal magnetic field was adjusted so that magnetization of 
the FL is perpendicular to the magnetization of the RL 
(6 = 90°). The magnitude of the RF input signal, con- 
nected to the MTJ through the capacitive lead of a bias 
tee, was fixed to -15 dBm. This resulted in a RF current 
(Irf) between 5 fiA and 25 /iA, depending on the sample 
resistance. Irf was calculated on the basis of the non- 
resonant background signal, using a model proposed in 
Ref. [8| . The bias voltage was fed through the inductive 
lead of the bias tee. y m i X was measured using a AC cou- 
pled lock-in amplifier, which was synchronized with the 
amplitude modulated signal from the RF generator. In 
this paper, positive bias voltage indicates electron trans- 
port from the bottom RL to the top FL. 



III. RESULTS AND DISCUSSION 

Table [^summarizes the TMR, the RA product and the 
static offset magnetic field (Hs) for three sets of MTJs 
with different MgO tunnel barrier thickness. The rep- 
resentative TMR vs. magnetic field loops are presented 
in Fig. □ The high TMR ratio of 170% for a 1.01 nm 
thick barrier and the exponential decrease in RA product 
with decreasing MgO thickness confirm good tunnel bar- 
rier quality [1J|. Similar TMR ratios and RA products 
were measured on full wafers using a current in-plane 
tunnelling (CIPT) technique before patterning [l2j |. The 
overall offset field (Hs) is shifted approximately 30 - 40 
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FIG. 1: TMR vs. magnetic field loops of samples S1-S3. 
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FIG. 2: ST-FMR spectra of samples SI (a), S2 (b) and S3 (c) 
measured with various magnetic field applied at an angle of 
f3 = 70° with respect to the easy magnetization axis. Only 
the RF signal (without DC bias voltage) was supplied to the 
MTJ. For sample S3 (c) two closely spaced peaks are visible. 
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FIG. 3: The dispersion relation of sample S3 measured with 
the magnetic field applied at an angle of /3 = 70° (a) and /3 
= 30° (b) with respect to the easy magnetization axis. The 
solid and dashed lines represent macrospin simulations of the 
FL and RL, respectively, (a) At an angle of /3 = 70°, the 
resonance frequency of two slightly separated FL modes in- 
crease with increasing magnetic field, (b) At an angle of /3 
= 30°, magnetization precessions in both the FL and RL are 
measured. 



Oe with respect to the wafer-level measurements due to 
dipolar magnetostatic stray-field coupling in the nanopil- 
lar junctions. For the MTJ with a 1.01 nm thick tunnel 
barrier, antiferromagnetic stray-field coupling dominates 
the interaction between FL and RL (Hs — -21.7 Oe). 
A reduction of the barrier thickness to 0.76 nm reverses 
the sign of the offset field (Hs — 47 Oe). In this case, 
the FL and RL couple ferromagnetically due to direct 
interactions across the thin MgO tunnel barrier. 



A. ST-FMR 

Typical ST-FMR signals (without DC bias voltage) for 
samples SI - S3 are presented in Fig. We note that 
a single symmetric peak is measured for sample S2 in 
a wide magnetic field range. For this sample, the cou- 
pling between FL and RL is negligible. Moreover, the 
monotonic increase of the resonance frequency with ap- 
plied magnetic field indicates that the FMR signal orig- 



inates from magnetization precession in the FL [jjj . A 
similar behavior is observed for sample SI, wherein the 
effective coupling between FL and RL is weakly antiferro- 
magnetic. However, for sample S3, which is characterized 
by strong ferromagnetic coupling between FL and RL, an 
additional peak is measured. The origin of this double 
resonance mode is not entirely clear. In previous publica- 
tions, it has been attributed to domain formation in the 
FL [lfl , higher-order spin wave excitations [lj| and mag- 
netization precession in other layers of spin- valve MTJs 
17]. To analyze the double resonance mode in sample S3 
in more detail, we performed macrospin simulations us- 
ing the model presented in Ref. [l8j]. This model, based 
on the Stoner-Wolfarth approach, assumes coherent ro- 
tation of the FL and RL magnetization. By minimizing 
the system energy we find the angle of the FL and RL 
magnetizations with respect to the easy axis and on this 
basis, we calculate the dispersion relation. The simulated 
dispersion relations that are obtained for (3 — 70° and for 
(3 = 30° are presented in Fig. [3] together with the mea- 
sured ST-FMR spectra. For /3 = 30°, the experimental 
and simulated FMR modes of the FL and RL are in good 
quantitative agreement. We note that the FMR signal of 
the RL is only measured when a large positive magnetic 
field is applied to the nanopillar junctions. The reso- 
nance frequency of the RL decreases with increasing field 
strength in this field range. The frequency of the double 
resonance peak in the spectra for j3 = 70° (Fig. [Ha)), 
on the other hand, increase with applied field strength. 
The experimental dispersion relations now closely match 
simulated curve. Based on this analysis, we attribute the 
double resonance mode to inhomogeneous magnetization 
precession in the FL rather than FMR in the RL or any 
other magnetic layer of the MTJ stack. 

To further elucidate the origin of double- mode FL spec- 
tra, we simulated the resonance characteristics of MTJ 
nanopillars using OOMMF software [lj| with an additional 
extension enabling calculations of TMR and STT effects 
[20| . In these micromagnetic simulations, elliptical mul- 
tilayer systems with a 2 nm thick FL, a 1 nm thick MgO 
tunnel barrier, a 2 nm thick high-anisotropy RL, antifcr- 
romagnetically coupled to a 2 nm thick exchange-biased 
pinned layer (PL), were used. The area of the junction 
was identical to the experimental structures. The inter- 
layer exchange coupling and anisotropy energies were ex- 
perimentally determined by magnetic and magnetotrans- 
port measurements. Variation of the ferromagnetic inter- 
layer exchange coupling from to 19 /LtJ/m 2 in the sim- 
ulations yielded results comparable to the experimental 
data. We note that dipolar coupling between the FL and 
RL is intrinsically calculated and taken into account in 
OOMMF. Thus depending on the strength of the interlayer 
exchange coupling (input parameter), the effective cou- 
pling between FL and RL varies from antiferromagnetic 
to ferromagnetic in accordance with the experimental re- 
sults on samples SI - S3. 

The dynamic simulations were conducted in the fol- 
lowing way: first, an external magnetic field was applied 
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FIG. 4: Simulated ST-FMR curves for weak (a) and strong 
(b) ferromagnetic interlayer exchange coupling. In OOMMF 
simulations, a voltage step was used to excite magnetization 
precession in the FL of a MTJ structure. The dimensions of 
the simulated and experimental junctions are identical. The 
existence of closely-spaced double-peak ST-FMR signal for 
strong coupling is independent of the anisotropy constant. 



at an angle with respect to the magnetic easy axis. Af- 
ter relaxation, a voltage step was applied to exert a STT 
on the FL. The voltage step amplitude was adjusted, so 
that the FL magnetization oscillations changed the MTJ 
resistance by a few Ohms. The used values correspond 
to an AC current of a few /iA, which closely mimic the 
experimental conditions and ensures that the magneti- 
zation oscillations are within the linear regime. Finally, 
the resonance spectra were obtained by Fourier transfor- 
mation of the time-derivative damped oscillation of the 
simulated tunneling magnetoresistance. 

Figure 2] presents the simulated ST-FMR spectra for 
two MTJ nanopillars that closely resemble experimental 
samples S2 and S3. The simulations confirm that the 
magnetization of the RL does not precess under these 
conditions (f3 = 70°) in the investigated frequency range. 
For a weak interlayer exchange coupling energy of J = 
6 /xJ/m 2 (sample S2), a single resonance peak is sim- 
ulated for different FL anisotropy energies - Fig. 2k - 
and different magnetic field strength (not shown), which 
fulfills the Kittel dispersion relation. For a larger ferro- 
magnetic coupling energy of J = 19 /jJ/m 2 (sample S3), 
an additional broad resonance peak was resolved in the 
simulations (Fig. SJb)), regardless of the FL magnetic 
anisotropy. This behavior is reminiscent to the experi- 



mental behavior of sample S3 with a 0.76 nm thin MgO 
tunnel barrier. The simulations thus confirm that the 
the double resonance mode originates from inhomoge- 
neous magnetization precession in the FL of the MTJ 
nanopillar stack due to strong interlayer exchange cou- 
pling between FL and RL. 



B. Torques and torkances 

In order to obtain the STT components, i.e., in-plane 
torque rii and perpendicular torque r_i_, from the ST- 
FMR measurements, we used the model presented in Ref. 
[l3[. Here, we assume a simplified formula for V m i x : 



V m 



l d 2 V 
4~dP J 
1 d 2 V 



UF 



(2a) 



H"/ sin 9 



2 dldO AeM s Vola 



RF 



(2b) 



where h is the reduced Planck's constant, 7 is the gyro- 
magnetic ratio, e is the electron charge, Vol is the volume 
of the FL, Ms is the saturation magnetization of the FL, 
a is the linewidth, £|| = 2(e/hsin6)(dV/ dI)d,T\\/ dV and 
£_l = 2(e/hsm,6)(dV/ dl)dr±/ dV are the magnitudes of 
the symmetric S(ui)=[l+(uj-ui m )' 2 /a 2 } -1 and asymmet- 
ric A((jj)=[(u)-u} m )/o~]S(uj) lorentzians components, and 
n ± = 1 N x M cS /uj m , N X =4:TT+{H Z -H a sm 2 (3)/M e s, where 
u> m is the resonant frequency, H z is the sum of the applied 
external magnetic field and the offset field acting on the 
precessing FL, H a is the in-plane anisotropy field of the 
FL and A-kM c s is the effective out-of-plane anisotropy of 
the FL. We neglected the terms (2c) - (2g) of Ref. [HI 
because in our case w 90°. 

Figure [5^ presents a comparison of the in-plane 
torkance in samples SI, S2, and S3. The absolute value 
of the in-plane torkance increases with decreasing barrier 
thickness and it only weakly depends on DC bias voltage. 
According to Slonczewski's free electron model for elas- 
tic tunneling in symmetric MTJs, the in-plane torkance 
is proportional to the differential conductance measured 
for parallel alignment of FL and RL [2l| : 



dV 



2p 



2e 1 



dl 
dV 



(3) 



By using Jullieres model to derive the spin polariza- 
tion of the tunneling current p at V = V, we found a 
good match between our experimental data and theoret- 
ical calculations based on Eqj3](Fig. [3a)). The absolute 
torque values in Fig. [SJb) were obtained by numerical 
integration of the data in Fig. [5ja). Obviously, the in- 
plane torque varies linearly with DC bias current and it 
is independent of MgO tunnel barrier thickness. These 
results are in good agreement with previously published 
experimental data in Refs [1, [j| [T3I I22II and calculations 
based on an ab initio approach [23j, |24J . 
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FIG. 5: Bias dependence of the in-plane torkance (a), in- 
plane torque (b), perpendicular torkance (c) and perpendic- 
ular torque (d) for MTJs with different MgO barrier thick- 
ness. The solid lines in (a) represent calculations based on 
Eq. [3] The torque values are numerically integrated from 
experimentally determined torkances. t± for sample S3 was 
compensated for an error originating from asymmetric ST- 
FMR resonances. 

Experimental data on the perpendicular torkance are 
summarized in Fig. [SJc). For samples SI and S2, the 
torkance decreases with DC bias voltage and dT±/dV — 
for zero DC bias voltage as predicted by theoretical cal- 
culations. However, a discrepancy is observed for sample 
S3. In this sample, strong ferromagnetic coupling be- 
tween the FL and RL of the MTJs results in asymmet- 
rical double resonance modes in the ST-FMR spectra. 
The fitting procedure based on Eq. [2] therefore intro- 
duces an error in the experimental torkance values for 
this sample. A good match with theoretical calculations 
is obtained when this artifact is compensated by subtrac- 
tion of a constant torkance value. Figure [SJd) illustrates 
that the absolute perpendicular torque varies quadrati- 
cally with DC bias current. Moreover, t± is similar for all 
samples. We note that different torque versus bias depen- 
dencies have been measured recently. Especially, it has 
been shown that the sha pe of r i (V) curves can change 
from quadratic to linear [2£|,|2q. However, such effects 



were only measured in asymmetric MTJs with different 
FL and RL electrodes. In our junctions, the composition 
and thickness of the CoFeB electrodes are the same 



IV. SUMMARY 

In summary, we have investigated MTJ nanopillars 
with varied MgO tunnel barrier thickness using the spin- 
torque diode effect. We measured a symmetric ST-FMR 
signal for samples with ^MgO > 0.9 nm. In this case, 
the coupling between FL and RL is weakly antiferro- 
magnetic. Contrary, double and closely-spaced resonance 
modes were obtained for MTJs with a 0.76 nm thick 
tunnel barrier. Macrospin and micromagnetic simula- 
tions indicate that the asymmetric double-peaks orig- 
inate from inhomogeneous magnetization precession in 
the FL caused by ferromagnetic coupling to the RL. The 
in-plane and perpendicular torques scale with DC bias 
current and they are independent of MgO tunnel barrier 
thickness. 
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